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I 1 INTRODUCTION 
• i— I . 

■ X-ray emission from late-type M dwarfs has been studied 
S_i ' extensively to investigate the structure and emission mech- 
anisms of stellar coronae. Coronal heating to X-ray emitting 
temperatures is attributed to either impulsive flares or qui- 
escent energy release in magnetic structures. The corona of 
these stars are thought to be similar to the sun but often 
with luminosities orders of magnitude higher. The high mag- 
netic activity of these flare stars is also seen in their optical 
spectra. Continuum enhancement and strong emission lines 
of the H Balmer series, Ca II and neutral He are often evi- 
dent (Montes et ah, 1999). 

An EXOSAT study (Pallavicini et al. 1990) showed 
flares have a wide range of energies and timescales. Most 
outbursts can be described as either implusive (decay time 
< 1 hr) or long decay flares (decay time > 1 hr) and have 
thermal X-ray spectra with temperatures similar to solar X- 
ray flares. The impulsive stellar flares have timescales simi- 
lar to the compact solar flares. The long duration flares have 
greater total energy and are more similar to two- ribbon flare 
events. 

From ROSAT observations, coronal emission from dMe 
stars has been shown to have two distinct spectral compo- 
nents, a low temperature component attributed to quiescent 



ABSTRACT 

We report the serendipitous discovery of a flare star observed with the ROSAT X- 
ray observatory. From optical spectra, which show strong and variable emission lines 
of the hydrogen Balmer scries and neutral helium, we classify this object as a M3.0Vc 
star, and estimate a distance of 52 pc from published photometry. Due to the star's 
close proximity (13.6' ) to the calibration source and RS CVn binary AR Lacertae, 
long term X-ray coverage is available in the ROSAT archive (^50 hours spanning 6.5 
years). Two large flare events occurred early in the mission (6-7/1990), and the end of 
a third flare was detected in 6/1996. One flare, observed with the Position Sensitive 
Proportional Counter (PSPC), had a peak luminosity La'=1-1x10 30 erg s _1 , an e- 
folding rise time of 2.2 hours and a decay time of 7 hours. This decay time is one of 
the longest detected on a dMe star, providing evidence for the possibility of additional 
heating during the decay phase. A large HRI flare (peak Lx=2.9xl0 30 erg s" 1 ) is 
also studied. The "background" X-ray emission is also variable - evidence for low-level 
flaring or microflaring. We find that >59% of the HRI counts and >68% of the PSPC 
counts are due to flares. At least 41% of the HRI exposure time and 47% of the PSPC 
are affected by detectable flare enhancement. 
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active regions and a variable, high temperature component 
due to compact flaring regions (Giampapa et al. 1996). 

Schmitt (1994) has shown conclusively the existence of 
long duration flares on M stars using the ROSAT all-sky sur- 
vey. In some flaring stars, the long decay can be attributed 
to continual heating of the flaring region (Schmitt & Fa- 
vata, 1999; Ottmann & Schmitt, 1996). New flare models 
have been developed (Reale et al., 1997) in which the ad- 
ditional heating determines the characteristics of the decay. 
Using this model, an analysis of long duration flares on AD 
Leo (Favata, Micela & Reale, 2000) and EV Lac (Favata et 
al., 2000), have shown the emitting regions to be compact 
with length scales less than the stellar radius and similar 
in size to solar flares, thereby providing evidence that long 
duration flares are produced in high pressure structures. 

During an analysis of observations with the ROSAT 
High Resolution Imager (HRI) of the RS CVn binary AR 
Lac, we noticed an X-ray source within the field of view to 
be highly variable and undetected in many fields. The first 
X-ray detection of this source was with the Einstein observa- 
tory (hence the catalogue name, 2E 2206.6+4517; Harris et 
al., 1993). Coaddition of six IPC observations spanning 26.6 
ksec yield a 4a detection with 115 net counts. The source is 
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undetected in a 1.5 ksec HRI observation. No flaring activity 
is evident during these observations. 

We analyzed thirty nine observations Irom the ROSAT 
public data archive (16 PSPC; 23 HRI) which included this 
source position within the field ol view. The total observing 
time was «50 hours. The observations span two flares with 
almost complete light curves, the tail end ol a third flare, 
and show variability in the low-level X-ray emission. The 
flare detected with the PSPC has a long decay time, possibly 
providing evidence lor significant continual heating during 
the flare decay. Thus, a lurther in-depth study ol flares and 
quiescent emission from 2E 2206.6+4517 could provide a 
uselul test ol current models. 



2 OBSERVATIONS AND DATA REDUCTION 

We detect the X-ray source 2E 2206.6+4517 in 16 obser- 
vations with the ROSAT Position Sensitive Proportional 
Counter (PSPC) and 23 observations with the HRI, within 
the instrument bandpass ol 0.1 - 2.4 keV. Extensive and 
continuous X-ray observations («17 hrs) were made with the 
PSPC between 18-22 June 1990 (during the ROSAT in-orbit 
calibration period), 30-31 December 1991 , and 29 May 1993 
and 02 June. The HRI calibration observations ol AR Lac 
began thirteen days alter the completion ol the PSPC ob- 
servations. Fourteen HRI pointings between 2-8 July 1990, 
include the source 2E 2206.6+4517. Eight additional HRI 
observations are available in the archive over the period ol 
June 1992 through November 1996 for a total observing time 
ol 33 hours. 

We measured count rates using the IRAF/PROS data 
analysis software and corrected for vignetting due to the 
large range ol off-axis angles (2' < 9 < 47' ). The observa- 
tions were subdivided into multiple time bins to achieve a 
higher temporal resolution while preserving a minimal 2a 
detection for each bin. An annular region was centered on 
the source to correct for the background count rate for most 
cases. For detections near the edge ol the field, a nearby cir- 
cular background region was chosen. A log ol the ROSAT 
observations is given in Table [l] which includes the exposure 
time and off- axis angle. 

Spectral fitting ol the PSPC data was done with the 
XSPEC software package. Source and background counts 
were extracted using the XSELECT task Irom the FTOOLS 
package, and the ancillary response files were constructed 
with PCARF to account for off-axis vignetting. We ignored 
the lowest 11 spectral energy bins to avoid scattered solar 
EUV contamination. The highest 56 spectral channels were 
also omitted due to poor statistics, a result ol a marked 
decrease in the effective area ol the instrument at these en- 
ergies. The energy bin distribution oversamples the intrinsic 
spectral resolution ol the PSPC, so we grouped the bins by 
a lactor ol 5 to improve the statistics. 

Multiple optical spectra were taken by Perry Berlind 
with the Tillinghast 60" telescope and FAST spectrograph 
(Fabricant et al., 1998) at the Fred Lawrence Whipple Ob- 
servatory on Mount Hopkins. A slit width ol 3" , a 300 
lines/mm grating, and the Loral CCD with 15/^m pixels 
provided a resolution of ~ 5 A. A 3.5 minute exposure was 
taken on UT 16 May 1998 and two ten minute exposures 
were acquired on UT 30 May 1998 and UT 24 June 1998. 
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Notes: rp=PSPCC; rp^PSPCB; rh=HRI 



An observation ol Feige 34 was used lor extinction correction 
and flux calibration. Standard bias subtraction, flat-fielding, 
the extraction ol one-dimensional spectra, and wavelength 
calibration were performed in IRAF. 



3 SOURCE IDENTIFICATION AND OPTICAL 
SPECTRA 

After obtaining an X-ray source position from a nearly on- 
axis HRI observation, we identified two candidate optical 
counterparts within 10" on the Digitized Sky Survey red 
plates. 

Two short observations were adequate to classify can- 
didate 1 as an active Me star based on its Balmer emission, 
and candidate 2 as a normal G dwarf, with no evidence for a 
composite spectrum. We thus identify the Me star as the op- 
tical counterpart to 2E 2206.6+4517. Coordinates from the 
USNO-A2.0 catalog (Monet et al., 1998) are a =22:08:37.5, 
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Figure 1. Digital Sky Survey plate showing the two possible 
optical counterparts to the X-ray flare. The star labeled dMc is 
the likely counterpart. The field of view is 5' x5' (North is up, 
East is left). The thin structures to the SE are plate artifacts. 
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Figure 2. Spectrum of the optical counterpart of 2E 
2206.6+4517 observed on 30 May 1998. The strong Balmer emis- 
sion lines are indicative of a dMe star in an active flare state. 



S = +45:31:27.9 (J2000), and a 5' finder chart is provided 
in Figure [jj 

To facilitate an accurate spectral classification, the 
emission-line star was observed again on UT 30 May 1998 
and UT 24 June 1998 for 10 min each. Both spectra showed 
strong emission lines, typical of a dMe flare star. 

The first long spectrum (Figure ^) showed strong emis- 
sion lines, typical of a dMe flare star in an active flaring 
state. The second observation (Figure ^) showed a marked 
decrease in the strength of the emission lines, indicating that 
the flare star was either in or close to quiescence. Table 
|i] lists the measured emission line equivalent widths. The 
strength of chromospheric heating due to the flare is evident 
by the increase in the hydrogen line emission by a factor of 
1.6 (H a)- 4.1 (H 7). 

The optical magnitudes from the USNO catalogue are 
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Figure 3. A second optical spectrum of 2E 2206.6+4517 was 
obtained on 24 June 1998. The decrease in the emission line 
strengths, as compared to the spectrum obtained one month ear- 
lier, could be attributed to a decrease or lack of flare activity. 



Table 2. Optical emission line equivalent widths 



Date 


30 May 1998 


24 June 1998 
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Note: a continuum difficult to define 



B=16.5 and R=14.2 for this source, colors consistent with 
a late-type star. However, these colors may be affected by 
the star's strong emission lines. To derive an accurate spec- 
tral classification, we first removed the emission lines from 
the spectra by interpolating to the adjacent 'continuum'. 
We then cross-correlated the resulting spectrum (using the 
IRAF fxcorr task) with a sample of digital spectra of M 
dwarfs from the Gliese catalog (Henry, Kirkpatrick, & Si- 
mons, 1994). By far the best correlation, and a very close 
match was with the M3.0V star Gliese 251 (cross-correlation 
peak height 0.99). 

For an M3.0V star, we find M v = 11.7 + 0.6 from equa- 
tion 1 of Henry et al. (1994; note that this error is the empir- 
ical rms value, more representative of the actual error than 
lcr). Bessell (1991) tabulates V — R = 1.1 and B — V — 1.55 
for an M3.0V star. From this we derive Mr — 10.6 for com- 
parison to the USNO red magnitude, to obtain a distance of 
52 pc. 



4 ROSAT X-RAY OBSERVATIONS 
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4.1 Flare Activity-Light curve 

The PSPC light curve (Figure shows prominent X-ray 
flare activity on 19 June 1990. Using the parameters for 
the best fit spectral model (section 4.3), the peak flare lu- 
minosity is Lx = l.lxl0 30 erg s _1 , a 24 fold increase over 
the mean quiescent luminosity Lx=4.6xl0 28 erg s _1 (Sec- 
tion 4.2). The outburst can be characterized by an e- folding 
rise time T r —2.2 hours and decay time Td—7 hours. The to- 
tal rise and decay time for the flare is At r i se — 6 hours and 
Atdecay— 30 hours as measured from the quiescent to peak 
count rate. The tail end of the flare significantly departs 
from the exponential decay of the flare (Figure W). 

This flare is evidently a long duration event. Most long 
decay flares on dMe stars have Td ~ 1 hour (Pallavicini et al., 
1990). Recently, long duration flares have been detected on 
EV Lac (T d =10.5 hours; Schmitt 1994) and AD Leo (r d =2.2 
hours; Favata, Micela & Reale, 2000). Continual heating of 
the flaring region during the decay has been proposed to 
explain such long decaying events. Due to the long decay 
time, the total energy Etot=3.5x 10 34 erg (Table 3) re- 
leased is similar to the flare seen on EV Lac (Etot=9x 10 33 
erg; Schmitt 1994) and large as compared to other dMe flare 
stars (^ 3 x 10 30 - 1 x 10 34 ; Pallavicini et al., 1990). These 
long decay flares on dMe stars are, however, 2-3 orders 
of magnitude less energetic than giant X-ray flares on RS 
CVn stars such as Algol (Td— 8.4 hours; Etot=7x10 36 erg; 
Ottmann & Schmitt, 1996), and CF Tucanae (1^=22 hours; 
Etot=14x10 37 erg; Kiirster & Schmitt, 1996). 

On 04 July 1990 at 21:11:35 UT, a flare with a peak 
luminosity, about the same magnitude as the PSPC out- 
burst (Figure was detected with the HRI. Since the 
HRI has extremely limited spectral resolution, we input the 
model fit to the PSPC flare and the HRI count rate into 
the IRAF/PROS task hxflux to convert counts to luminos- 
ity. The peak flare luminosity was Lx=2.9x 10 30 erg s _1 
a factor of 54 larger than the mean quiescent luminosity 
Lx=5.4xl0 28 erg s _1 . The outburst can be characterized 
by an e- folding rise time r r ~15 minutes and decay time 
Td— 1.2 hours. The total rise and decay time for the flare is 
Atrise— 40 minutes and Atdecay—4.6 hours. The total energy 
released during the flare is Etot=L6x 10 34 erg (Table 3). 

We also discovered the tail end of an additional flare ob- 
served on 03 June 1996 with a two-fold rise above the back- 
ground level. Due to our incomplete light curve, we can only 
provide an approximation to the decay time ( ~9 hours). 
The X-ray characteristics of these flares are listed in Table 



I 



4.2 



Count Rate Distribution Analysis: Quiescent 
Level and Flaring Fraction 



The "quiescent" X-ray emission outside the large flares is 
clearly itself variable. Fig. 4 shows a factor of 3 change in 
the level outside of flares, and three long HRI observations 
between 12/95 and 7/96 show average of Lx = 15.4, 6.6, 
and 5.0 x 10 28 erg s" 1 . 

To investigate this further, we performed a statistical 
analysis of the complete X-ray light curve to determine the 
minimum fraction of counts due to flares, and the minimum 
fraction of time spent in a flaring state. We assumed a truly 
non-variable, quiescent background rate exists that can be 
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Figure 4. X-ray flare of 2E 2206.6+4517 observed with the 
PSPC. The flare begins on 19 June 1990 at 12:45:22 UT and lasts 
for 1-1.5 days. The error bars are ±ltr based on count statistics. 
The x error bars show the time spanned for each measurement. 
The dashed line is an exponential function characterized by the 
e-folding time, Td—7 hours. The soft and hard bands for the hard- 
ness ratio are defined as 0.07-0.42 keV and 0.42-2.48 keV. 
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Figure 5. X-ray flare of 2E 2206.6+4517 observed with the HRI 
on July 4, 1990. The flare was first detected at 21:11:35 UT and 
lasted for 4.6 hours. The lower plot focuses at the outburst in 
detail. The dashed line is an exponential function characterized 
by the e-folding time of the decay t,2— 1.2 hours. 



described by a Poisson distribution. We then determined 
the fraction of counts due to flares by integrating the counts 
left unexplained by a least squares fit of a Poisson func- 
tion to the low end of the count rate distribution (Saar & 
Bookbinder 1998) . We averaged the results of fits using sev- 
eral reasonable count rate binning sizes. This analysis yields 
the minimum fraction of flare counts, since flares below the 
instrument sensitivity, and possible quiescent level changes 
due to rotation and evolution of magnetic regions are all in- 
cluded in the derived quiescent level. We find that (at least) 
68 ± 2% of the PSPC counts and 59 ± 2% of the HRI 
counts are due to flares. Upper limits to the quiescent flux 
level and the fraction of time that the star is quiescent are 
also derived from the analysis. We find that the (Poisson) 
mean quiescent count level is 0.022 ± 0.001 counts s _1 for 
the PSPC {L x ~ 4.6 ± 0.2 x 10 28 erg s _1 ) and 0.0056 ± 
0.0006 counts s^ 1 for the HRI (Lx ~ 5.4 ± 0.6 x 10 28 erg 
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s 1 ). The minimum fraction of time with a detectable flare 
contribution to the observed flux is 47 ± 5% (PSPC) and 
41 ± 8% (HRI). 

All of the PSPC data was taken between 1990 and 1993, 
while the HRI data include a significant fraction from 1995 
and 1996, with an exposure time weighted time difference of 
~1200 d. Thus it is possible that the small (~15%) difference 
in the PSPC and HRI quiescent fluxes may be explained by 
time evolution of the quiet emission, due to e.g., long-term 
evolution in the numbers of active regions. Unfortunately, 
there are relatively few quiescent counts in the 1995-96 HRI 
data, making a direct test of time variation inconclusive. 
Further data over a longer timescale would help decide the 
level of the star's quiescent coronal variability. 

We find that the quiescent and bolometric luminosity 
agree with the linear correlation between these quantities, 
as noted by previous studies of flare stars (Pallavicini et 
al., 1990; Agrawal et al., 1986). We estimate a bolometric 
luminosity of 6.74xl0 31 erg s -1 , using M v = 11.7 and the 
bolometric correction of Pettersen (1983). 

4.3 Spectral Analysis 

We extracted spectra from the brightest quiescent and flare 
PSPC images. During the longest pointing of this field, the 
source 2E 2206.6+4517 was in quiescence (rpl00588) and 
16.5' off-axis. The data set corresponding to a flare event 
with the highest number of source counts was rpll0591, for 
which the source was nearly on-axis (4.4' ). 

For the quiescent X-ray emission, we found that no 
single-component thermal plasma (Raymond and Smith 
1978) model could adequately fit the data. Unfortunately, 
the signal to noise of our spectrum was not sufficient to 
uniquely constrain a two-component model. However, we 
found that models that did not contain a Raymond-Smith 
component of kT w 1.0 keV could be rejected (i.e. xt » !)■ 
These models significantly underpredict emission from ap- 
proximately 0.85-1.0 keV, corresponding to the Fe L-shell 
blend. However, with an appropriate "high" temperature 
thermal plasma model, the remaining low temperature emis- 
sion has xt statistics much less than one, indicating that the 
models are not well constrained. Despite the poorly con- 
strained two component model, the spectral fit suggests the 
existence of a coronal plasma of kT ~ 1 keV. A compari- 
son of the normalizations of the two models shows that the 
"hot" and "cool" components contribute nearly equal flux 
in the PSPC band (Table g). 

Because of the fewer counts in the flare observation, a 
larger binning factor was necessary to provide a significant 
S/N for spectral fitting. However, we could not fit the result- 
ing spectrum with any one or two component model. The 
heavy binning broadens the effect of uncertain calibration 
features, particularly the PSPC window carbon edge at 0.4 
keV. To compare the quiescent and flare spectra, we must 
mitigate this effect. The peaks and troughs in the spectra 
most affected by binning all occur below about 1 keV. Using 
only the photons above 1 keV has the added advantage that 
the "hot" component dominates in this regime, so we may 
use a single temperature model. We fit both spectra in the 
energy range 1.0-1.8 keV with single temperature Raymond- 
Smith model, using the value of Ah from the full quiescent 
spectrum. 



A direct comparison of the flare and quiescent normal- 
ization shows an increase by nearly a factor of 20 in emission 
measure. In addition, though the temperatures are not very 
well constrained, their la errors just barely overlap, indi- 
cating that the increase in X-ray emission was likely accom- 
panied by an increase in the plasma temperature. This is 
consistent with behavior seen during flares in late-type ac- 
tive stars ( e.g., Giampapa et al., 1996, Singh et al., 1999). 

The small number of counts does not allow a determi- 
nation of a temperature variation during the flare decay. To 
investigate the spectral evolution of the flare, we calculated 
hardness ratios (Figure 4). The count distribution appears 
to harden by about 40 — 50% during the onset of the flare. 
A decrease of the hardness ratio during the decay is evident 
without certainty due to limited count statistics. 

Singh et al. (1999) found that Fe abundances for ac- 
tive dwarf stars with log (Lx/Lboi) > —3.7 are strongly 
subsolar, whereas the Fe abundances in the less active stars 
are within a factor of 2 of the solar value. Since the quies- 
cent luminosity ratio for 2E 2206.6+4517 hovers near this 
value, low abundances may also affect our X-ray spectral 
fits. On the other hand, flare activity has been observed 
to increase the apparent elemental abundances (Ottman & 
Schmitt 1996). 



5 CONCLUSION 

In this paper, we have identified a new flare star 2E 
2206.6+4517 using extensive ROSAT X-ray observations 
and optical spectra. The object has been classifed as a 
M3.0Ve star. Three large X-ray flare events are seen with 
light curves characteristic of outbursts detected in other flare 
stars. Two flares with almost complete X-ray coverage show 
a variation in strength and timescales. A flare detected with 
the PSPC had a peak luminosity Lx = l.lxl0 30 erg s _1 , an 
e-folding rise time of «2.2 hours and a decay time of f»7 
hours. We interpret this long decay time (one of the longest 
ever observed) as possible evidence for continual heating as 
similar to recent observations and analysis of AD Leo and 
EV Lac (Favata et al., 2000; Schmitt, 1994). An observation 
with the HRI detected a flare with a higher peak luminos- 
ity of Lx=2.9xl0 30 erg s , an e-folding rise time of ~15 
minutes, and a decay time of «1.2 hours. 

We used a statistical analysis of the X-ray light curves to 
measure the quiescent X-ray level, the (minimum) fraction of 
counts in flares and (minimum) fraction of time that the star 
exhibits flare or microflare activity. The PSPC data showed 
a quiescent luminosity of Lx = 4.6 + 0.2 x 10 28 erg s _1 , with 
>68±2% of counts coming from flares, and significant flare 
contribution to at least 47+5% of the time observed. The 
HRI data, taken on average almost three years later, show a 
quiescent luminosity of Lx = 5.4 ± 0.4 x 10 28 erg s^ 1 , with 
>59±2% of counts coming from flares, and significant flare 
contribution to at least 41+8% of the time observed. 

We obtained two optical spectra that show strong and 
variable emission lines of the hydrogen Balmer series and 
neutral helium which is further evidence for flare activity. 

With the large amount of X-ray coverage of the source 
2E 2206.6+4517 in the ROSAT data archive, this object 
provides an excellent opportunity to facilitate in the under- 
standing of the physical properties of these flare stars in- 
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Table 3. X-ray Flare Characteristics 

Instrument Date Peak(UT) At r i ae r r 

PSPC 19 June 1990 12:45:22 ~ 6 hrs ~ 2.2 hrs 

HRI 4 July 1990 21:11:35 ~ 40 min ~ 15 min 

HPJ 3 June 1996 



At 



decay 

30 hrs 
4.6 hrs 

9 hrs 



Td 



7 hrs 
1.2 hrs 



Lx(peak) 

l.lxlO 30 erg s" 
2.9 xlO 30 erg s" 



E T QT 

~ 3.5 X 10 34 erg 
~ 1.6 x 10 34 erg 



Table 4. X-ray Spectral Model Fits 



Model N H kT(l) Norm(l) kT(2) Norm(2) \l 

(cm" 2 ) (keV) (kcV) 



Quiescence RS + RS 5.0xl0 19 0.16" 4.8xl0" 5 0.97° 6.8xl0" 5 0.66 

Quiescence RS 5.0xl0 19 0.991 

lare RS 5.0xl0 19 1.321 

Note: (a) The one sigma error is not bounded 



Quiescence RS 5.0xl0 19 0.99l^ 6.8xl0" 5 0.33 

Flare RS 5.0xl0 19 1.32+g-| 1.3xl0 -3 0.66 



eluding the flare, quiescent and long term activity. Further 
measurements of this object may be afforded by calibration 
observations of AR Lac by Chandra or XMM-Newton. 
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